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Abstract 

The mean and the variance of the logarithm of the conductance {Ing) in the locahzed regime in 
the one-dimensional Anderson model are calculated analytically for weak disorder, starting from 
the recursion relations for the complex reflection- and transmission amplitudes. The exact recursion 
relation for the reflection amplitudes is approximated by improved Born approximation forms which 
ensure that averaged reflection coefficients tend asymptotically to unity in the localized regime, for 
chain lengths L = Na — > cxd. In contrast the familiar Born approximation of perturbation theory 
would not be adapted for the localized regime since it constrains the reflection coefficient to be less 
than one. The proper behaviour of the reflection coefficient (and of other related reflection param- 
eters) is responsible for various anomalies in the cumulants of In in particular for the well-known 
band center anomaly of the localization length. While a simple improved Born approximation is 
sufficient for studying cumulants at a generic band energy, we find that a generalized improved 
Born approximation is necessary to account satisfactorily for numerical results for the band center 
anomaly in the mean of \n.g. For the variance of In g at the band center, we reveal the existence 
of a weak anomalous quadratic term proportional to L^, besides the previously found anomaly in 
the linear term. At a generic band energy the variance of In (7 is found to be linear in L and is 
given by twice the mean, up to higher order corrections which are calculated. We also exhibit the 
L =independent offset terms in the variance, which strongly depend on reflection anomalies. 
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I. INTRODUCTION 

The advent of the scahng theory of locahzation in dimensional disordered systems |l| 
and more detailed developments of it in ID and for quasi ID systems [4| has inaugurated 
a golden age for mesoscopic physics, particularly the study of transport phenomena. 

The fundamental hypothesis in the scaling theory l| is that the scaling of the logarithm 
of a typical conductance g as a function of a characteristic size L of the system is described 
asymptotically for large L by a universal function, (3{\ng), of a single parameter (SPS), 
namely Ing itself, such that din g/ din L = j3{lng). The function j3 which may generally 
depend on dimensionality is independent of L and of microscopic parameters in the system. 
We recall that in the studies of scaling in ID systems the parameter ln(l + p), with 

p = i the resistance, was identified as the convenient scaling variable both in the localized 
regime (p >> 1) where it is self-averaging and in the low resistance (p << 1) (quasi-metallic) 
regime. This variable reduces to — Ing above for p >> 1 (localized regime) and, thanks to 
the Landauer formula, P = j£ (with and the reflection and transmission coefficients of 
the system, respectively), it coincides with — Inti — — In (7. 

In {2!, 3] it was argued that the scaling theory of Abrahams et al. [ll had to be interpreted 
in terms of the scaling of the distribution Pg{g) of the random conductance of the system. 
SPS means that Pg{g) is fully determined by a single parameter such as e.g. the mean 
logarithm, < Ing >, which is itself defined by a scaling equation of the above form with In^f 
replaced by < Ingf >. 

After years of debates the question of the validity of SPS in the theory of Abrahams et 

......-—.-.-".KaBBBa ,...„,..,.. 

justification of the SPS hypothesis in the analyses [2|, |3| rests on a random phase approxi- 
mation (RPA) which assumes that the phases of the amplitude refiection- and transmission 
coefficients Rl and Ti (with = \RL\'^,tL = \Tl\'^) are uniformly distributed over (0,27r) 
in the localization domain, i.e. for length scales L much larger than the localization length 
^. Despite the existence of strong evidence, both numerical [l^ and analytical H, 12| for 
uniform phase distributions for L >> ,^ in the ID Anderson model [l^ the SPS controversy 
is not resolved, essentially because phases and conductance are not independent random 
variables. 

Doubts about the validity of results based on RPA have recently led Deych and collab- 
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orators to reconsider the scaling problem for the exactly soluble Lloyd model [5j as well as 
to present simulation results for^he conductance distribution in the ID Anderson model in 
the region of fluctuation states [8']. On the other hand, Schomerus and Titov [gI, Q] (see also 
Roberts [14]) have discussed simulation results for the first four cumulants of Ing for the 
Anderson model for weak disorder, both at a generic band energy and at special energies 
(band center and band edges) where ordinary perturbation theory fails. They also discussed 
the cumulants using^ Fokker-Planck approach for the joint distribution of Ing and of the 
transmission phase p, Q]. Their results support the validity of the lognormal SPS form of 
the conductance distribution, at a generic band energy (unlike results of Roberts|l4|), while 
showing deviations from SPS at the special energies. 

Detailed analytical studies of conductance cumulants based on properties of symmetric 
groups defined from generalized transfer matrices, and on analytic continuation procedures, 
have been published earlier by Slevin and Pendry ISj and by Roberts [14]. However, the 
validation of these approaches rests on support from numerical simulations and studies of 
limiting cases [1 61] . 

In this paper we adopt a new direct approach, circumventing RPA, for studying the 
conductance distribution (the distribution of — In^f = — Int^.) analytically in the localized 
regime, for weak disorder. An essential ingredient of our analysis is the identification of a 
general type of anomalous (non-perturbative) effects in various complex reflection amplitude 
moments and reflection coefficient moments of a finite chain in the localized regime, both for 
a generic band energy and at the band center. These refiection anomalies strongly infiuence 
the logarithmic conductance cumulants. In particular, the anomaly of the second moment 
of complex refiection amplitudes at the band center is found to be responsible of the well- 
known Kappus-Wegner anomalies in the localization length 17, Q, 19] and in the variance 
of the logarithm, conductance y. 

In our analysis we choose to carry out the perturbation theory to fourth order in the 
disorder (i.e. the random site energies in the Anderson model). Now, by definition we have, 
for large L, {—Ing) = (— Int^) = 2L/^ (whith (...) denoting averaging over disorder), 
where the inverse localization length 1/^ is proportional to the variance of the site energies, 
e"^ = (el) {{en) = 0), for weak disorder |20]. It follows therefore that in the localized regime 
{L » ^) we have L e'^ » 1. Hence it follows that the n-th moment if — Int^, is of 
order {L e^)". This shows that the fourth order perturbation theory may yield a correct 
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description of the mean and the variance of —\ng only. While the study of the third- and 
fourth cumulants, respectively at sixth and eighth orders, is thus left for the future, our 
results for the first two cumulants at a generic energy already rule out single-parameter 
scaling if effects beyond leading order {Le'^) are retained. 

In Sect. II. A we present the analysis leading to exact formal expressions for the first 
and second moments of —Ing in terms of cumulated reflection amplitudes- and reflection 
coefficients moments. In II. B we define successively the improved Born approximation and 
a generalized improved Born approximation form of the exact recursion relation relating the 
random reflection amplitudes of samples of length n a and (n — l)a, respectively. These ap- 
proximate recursion relations are solved exactly to obtain explicit expressions for reflection 
amplitudes- and reflection coefficients moments. The generalized improved Born approxima- 
tion turns out to be vital for studying the band center anomaly in the localization length, in 
particular for achieving good ag reement with the numerical results of Kappus and Wegner 



and others 



17 



18 



191] . In Sect. Ill we discuss our detailed analytical results for 



the localization length and for the variance, both at a generic energy and at the band center. 
Some final remarks follow in Sect. IV. 



II. SCALING IN THE ANDERSON MODEL FOR WEAK DISORDER 
A. The mean and the variance of — Inti 

The Schrodinger equation for a chain of disordered sites 1 < m < of spacing 
a = l{L = N) is 

(fn+l + + £n^n = Eifn , (l) 

where the site energies in units of a constant hopping rate are mutually independent vari- 
ables, uniformly distributed between — f and -f (< e^+^ >= 0, < >= W'^/12, {e^J = 
W^/SO, etc). The disordered chain is connected as usual to semi- infinite non- disordered 
chains {em = 0) at both ends, with sites m < 1, and m > N, respectively. 

The distribution of the transmission coefficient t^r = I^tvP (conductance) for an electron 
incident from the right with wavenumber —k (energy E = 2 cos k) may be obtained by 
solving the general recursion relations which connect the complex transmission- (reflection-) 
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amplitudes Tn{Rn) of a chain of n sites with the corresponding amplitudes for a chain with 
one less disordered site, of length n — 1. These relations derived in|l2(] are, respectively, 

T = (2) 

" l-zz.„(l + e2^'=i?„_i) ' 



with 



The fundamental unitarity property. 



" 2 sin /c 



(4) 



|i?„P + |T„|2 = l , (5) 

follows quite generally from (j2H3]) by expressing using ([2]), and |-RnP, using ([3]), after 

rewriting the latter as i?„ = -1 + (1 + e^*''i?„_i)[l - iz/„(l + e^*''i?„_i)]~^ 
From ([2]) we obtain, with the boundary conditions Tq = 1 and -Rq = 0, 



n=l ^ 



N 

2ik ; 



2Z/„(1 + e'^'Rn-i] + c.c. , (6) 



which, together with ([3]), is our starting point for studying the probability distribution of 
— In t„ ~ — In via the calculation of its moments and the corresponding cumulants. We 
choose to restrict our analysis of the moments rrij = ((— IntAr)-^), j = 1,2,..., to effects up 
to 4th order in the random site energies and so we expand ([6]) in the form 



- In = 5^ 5^ ^^Hz/„(l + e2^^i?„_i)r + c.c. . (7) 

n=l p=l \ P ' 

As discussed in Sect. I systematic expansion to 4th order permits the explicit study of the 
first two moments only, in other words the determination of the inverse localization length 

1 fn\ 

- = —,N^oo , (8) 
^ 2iV' ' ^ ' 

and of the variance, 

var (— Intjv) = ^2 ~ N ^ oo . (9) 

After some calculations, using ([7]), the moments to 4th order in the explicitated z/„ are 
reduced to the following expressions (with z/p =< >): 



z/4 

mi = (z/2 )A^ + z/2 

2 



,2jA; p(l) ^ 



where 



ms = -u^{e^'^B!i^ + c.c.) + A^z/^ + N{N - l){u^f 



(10) 



(11) 



R^^ = Y,<Rl_,> , (12) 



n=l 



TV 



n=l 



The factorization of averages in fllOtilll) resuhs from the fact that the random amphtudes 
Rn-i and Rn-i are hnear functionnals depending on z/i, . . . z/„_i but not on Un, as shown 
by iterating ([3]), with Rq = 0. We have used the fact that the odd moments of Un are 
zero, which shows e.g. that < Um^^Rm-iRn-i >= 0,m n,p = 1,2. Also, in (11 OH 11 1) we 
have systematically dropped all terms of order higher than 4, in particular terms of the 
form < z/^i/^(_Rm-i-Rj^_i)^ >,P = l,2,m 7^ n, in ( ITT]) which are of orders 4 + 2p, unlike 
the corresponding terms with m = n which lead to lower order anomalous effects, as shown 
below. 

The Eqs. (llOtllip reduce the study of the mean and of the variance of — IntTv to the 
calculation of the quantities R^]^\ R^^^ and Qn\Qn\ which we refer to as fictitious cumu- 
lated reflection amplitudes- and cumulated reflection coefflcients moments (sums over chains 
having lengths equal to rational fractions of Na), respectively. These sums are dominated 
for — s> 00 by terms linear in A^, reflecting the fact that the amplitudes Rn for n — 00 
are described by an invariant (stationary) distribution 2l|, |2^ which is independent of the 
initial site where the iteration of ([3]) was started. The above cumulated moments will be 
studied in Sect. Ill, using the improved Born approximations discussed below. 



B. Improved Born approximations 

The Born approximation for the random reflection coefficient |-RnP for weak disorder is 
obtained by assuming i?„ to be typically proportional to z/„ and approximating ([3]) by the 
linear recursion relation Rjsi = e^''^Rn-i + whose solution 
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i?n = ^ J] e2*'=("-'"V™ , (14) 



m=l 



yields 



(|i?.r)=ni;2 = ^ , (15) 

using the familiar perturbation expression for the Anderson localization length|20|]. Since 
the absolute limit of (I-RtvP) for a chain of length is unity, it follows from f|T5|) that the 
Born approximation is not suited for discussing the strong localization (localized) regime, 
where Nv"^ » 1 oi N » The same conclusion also follows when using f|T^ to calculate 
the second moment of the reflection coefficient which enters in the definition of Q\f . In this 



(n-l)(z/2)2 + zy4 



which is also meaningless outside 



case, we obtain from f[T^ . (|-Rn|^) = n 
the perturbative domain {nv'^Y << 1- 

In contrast, the study of the localized regime is possible if one uses the improved (first) 
Born approximation where ([3]) is approximated, for weak disorder, by 

Rn = ''^'f-^. ^ + 0{ul) . (16) 

For example, by iterating the recursion relation for the averaged refiection coefficient, 
obtained from ( |T6l) (using the fact that Rn-i, is independent of i/„) and summing the resulting 
geometric series one finds 

{\Rn\') = l~{l-a,r,a, = {ul{l + ul)-') , (17) 

which has the desired limiting value of 1 for n ai » 1, while reducing, to leading order, to 
the perturbation result ( ITSjl in the opposite limit, n ai << 1. It follows that (|T6l) represents 
the simplest approximation of Eq. which permits a meaningful study of non-perturbative 
effects in the Anderson model in the localized regime. 

On the other hand, from the above discussion it is clear that the expression for the 
second moment of i?„ obtained from f|T^ . namely (R^) = —e^'^^^'nv'^ is also invalid in the 
localized regime. In contrast, by solving the recursion relation for (-R^) obtained from (fT6|) 
and summing the corresponding geometric series we get, for any 
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r 

(^n) = ^i^^ + ^(-')' = ^ = 0' 1' -'^"^ » 1 ' (18) 

neglecting exponentially small terms proportional to exp{n Inci). Here and in the following 
(9(z/^) refers to contributions from terms of pth order in the random site energies. Eq. (ITSll 
shows that while to leading order (-R^) is proportional to z/^ for 7^ (as in perturbation 
theory), it is strongly enhanced by non-perturbative effects, leading to (Rl) = — | + O^u'^) 
at E = 0. Using this value one obtains from ^ and ^ (with = 0{u^)) | = ^, 

which corresponds to a reduction of 33% of the result ^ = obtained by perturbation 



theory 20l|. This modification of the inverse localization length at the band center has the 



same origin as the well-known Kappus-Wegner|l7l| anomaly. But, clearly, its magnitude 
calculated within the improved first Born approximation is much too large since the effect 
obtained by Kappus and Wegner and by others 0, Q, Q, Il9| . using various sophisticated 



approaches, range between 8 and 9% of the perturbation result. This leads us to suggest a 
generalised improved Born approximation of Eq. ([3]) for dealing specifically with the Kappus- 
Wegner anomaly, which affects both mi, and 1712 at the band center, as shown by ( ITOl) and 
(ITTll . We refer to this more accurate procedure (for the band center) as the generalized 
improved Born approximation. In this approximation, besides the terms of the improved 
first Born approximation, we retain the term iVnC^'^^Rn-i in the numerator of ([3]) as well 
as a further term proportional to Rn-i obtained by expanding the denominator around the 
improved form (1 — ivn)^^- This yields the approximate recursion relation 



IV, 



Rn — QuRn-l + fn , fn — , Qn — JZ : , (19) 

1 - lUn (1 - lUnY 

where we have ignored a third order term proportional to Rn-i- In Sect. 3.2 we return to 
a crude estimate of the effects of this non-linear term to show that it does not affect the 
inverse localization length |, nor the variance of — Int^v, in the localized regime, to leading 
order in the disorder. 

III. LOCALIZATION LENGTH AND VARIANCE OF -Intjv 

We first analyse the moments of — Int^v at a generic band energy and then we discuss 
how the calculations are modified to account for the band center anomalies. 
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A. Generic band energy 



By averaging (fT6l) over disorder (using the fact that is independent of Rn-i) and 
iterating the resulting recursion relation for (Rn) in terms of a geometric series, which is 
readily summed, we get 

= u^u{N + u) + 0{u^),u = (e'^'^ - 1)-^ , (20) 

where we have ignored exponentially small terms proportional to e~^^^ for weak disorder. 
Proceeding in a similar way with the equation for (-R^) obtained from (fT6|l . we find 

R^^^ = u^v{N + v) + 0{u^),v = {e^''' -l)-^ . (21) 
From (pro]) and we then obtain 

m^ = {-\ntN)=(^^ + \[^{u^f-u^]^N-{u^Y{2\u\'' + \v\'') , (22) 

where the term linear in yields the familiar fourth order perturbation expression for the 
inverse localization length (1/^) at a generic energy jl^. I22I and the second term represents 
a new constant offset of mi. 

Next we obtain the asymptotic form of the cumulated reflection coefficient moments Q^'* 
and Q^"*. By summing f[T7l) over the disordered sites we find 

= AT - ATfli ^ 00 , (23) 

up to exponentially small terms. The term ^ is the leading deviation of Qtv"* from the 
unitarity limit (A^) in the localized regime, Nai ~ Nu"^ = ^ >> 1, in the improved Born 
approximation. 

In order to determine Q^^ we have to solve the two-point recursion relation for 
derived from (fT6ll . namely 



{\Rr,\') = 6o(|i?n-i|') + &2 + 46i(|i?„_iH - h^{e^'^{Rl_:, + c.c.) , (24) 
to 4th order non-vanishing terms. Here we have defined 

6, = (i.^^(l + i.^)~2),p = 0,l,2,... , (25) 
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and {\Rn-i\^) is given by ^ and {Rl_^) = (z/2(l - iz/„)-2)t;(l - e^^M^-i)) + 0{u^). As 
an example of the form of the solution obtained by iterating such a recursion equation we 
refer, for brievity's sake, to a similar equation which is solved in Scet. III.B. In the 4th 
order expression ffTTl) we require Q^^^ to negative orders in u up to zeroth order only. After 
performing successively the summation over sites in the solution of fl24l) and the further 
summation over sites in the definition ( fT3l) of Q)^ we obtain 



2bi + 62 



From (ED, dm), (I2TJI23D and (EMSE]) we then obtain, to order ^/^ 



^ ^ ^ ' ai 2(261 + 62)2 1-COS4A; ' ^ ^ 

which again involves a dominant term proportional to N and a higher constant offset 

term. The leading term 2v'^N in var (— Intjv) coincides with the result of various ear- 
lier theories , [14 . [15I 1 for a generic band energy. We recall that its proportionality to 
the dominant term of (— Intjy) in fl22]) ensures single-parameter scaling of the lognormal 
conductance distribution to lowest order, assuming that the higher cumulants are negligible. 
Our analysis reveals that this basic property is, in fact, a direct consequence of the unitar- 
ity limit of the reflection coefficient ([T71) (which defines Q^^ for asymptotic lengths in the 
localized regime). 

On the other hand, the cumulated refiection coefficients moments ( 123|) and ( 126|) are respon- 
sible for the existence of non-perturbative constant offset anomalies in the variance 0271) in 
the localized regime, for a generic energy. Finally, we observe that Eqs. (122|) and (127|) in- 
volving terms proportional to and additional constant offset terms conform to the ansatz 



of large- deviations statistics for cumulants % |23| . To zeroth order, the offset in the variance 



( |27|) reduces to the numerical constant —2. An analogous zeroth order offset constant of 
value — ^ has been obtained in|14i]. 



B. Band center 



At the band center the denominators in 0211) are singular which would require replacing 
them by their actual form e^'^^Ci — 1 in (fT8|) obtained from the improved Born approximation. 
However, as discussed in Sect. II. B we wish to further improve the calculation of at 
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the band center by using the generahzed improved Born approximation (fT9i) of the exact 
recursion relation Q. 

By squaring fll9p and averaging over the disorder we obtain the following relation for 
determining (-R^): 

{Rl) = e'''C{Rl_,) + 2e'''D{R^^r)-B , (28) 

where 

B = -{f^),C={gl),D={gM • (29) 

On the other hand, by averaging fll9l) . solving for (-R„) and performing the summation over 
sites in the solution we get (with / = {fn),g = {Qn)) 



g2«fc^ — 1 

By inserting fl30|) in fl28l) . the recursion relation for (-R^) takes the form 



^Rn) = ' , / . (30) 



{Rl)=e'^^C{Rl_,)-F{e'^'gY~' + F-B , (31) 

where 



1 — e^^'^g 

The exact solution of fl3T]) . with Rq = is given by 



F = ^Df-^-— . (32) 



n-l 



{Rl) = -BG''-^-F^G"-^{e^'^g)"' + {F-B)^G''-'^,G = e^'^C , (33) 



m=l m=2 

which, after performing the geometric sums, reduces to 



+ (^^-'^r'j^ - 1^ ■ (34) 

Finally, we evaluate R}^ defined by ( TT2l) and ( IMl) . ingoring exponentially small terms pro- 
portional to e^^^'^ (for weak disorder) in the localized regime. Specializing to the band 
center (fc = |) we have 
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R 



(2) 
N 



B 



Fg 



B-F 



+ 



N{B - F) 



(35) 



C{l-C)\ C + g l-C) (C-rg){l + g) l-C ' 

where the anomalous denominator 1 — C of order z/^ is responsible for the Kappus-Wegner 
correction in the localization length and a corresponding anomaly in var (— Int^r). By eval- 
uating the quantities entering in (l35l) and defined in ( |T9i) . ( l29l) and ( l32l) . for weak disorder, 
we obtain explicitly 



(2) 

AT 
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1 + 77^ + 3z/2 + (9(z/^ 



+ 



100 



(36) 



We recall that this quantity enters with opposite signs in (ITOl) and (ITT!) , respectively. 

The final expression of (— IntAr) at the band center is then obtained by substituting ([20 
and ( |36ll in (ITOl) . which yields 



+ 



100 



1 + 3z/2 + ^ + 0{u' 



(37) 



On the other hand, var (— Int^v) is obtained from ([9]) and (fTTj) by inserting (123|) and (!26ll . 
and (13611371) ■ This leads to 



1 Q - 

var(-lnt^) = ^(z.2)2iv2 + 



2.182Z/2 - 1.02(1/2)2 + o.Olz/4) 



N 



1 

Too 



2.01 + 6.061/2 + 8.08^ + 0{u^ 



2z;2^z^ (46^+362) _ ^3g^ 



ai 2 (261 + 62)' 



The inverse localization length at £" = obtained from ([3 



1 



1 

20 



9i72 + 7(^-2)2 _ 
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(39) 



may be compared to leading order with Thouless' perturbation expression l/,^o = z/2 /2j20| 



We thus find that the band center anomaly reduces the numerical coefficient of the z/2 term 
in 1/^ by 10% with respect to its value in Thouless' expression. For comparison, the earlier 
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n 



studies of the band center anomaly in the inverse locahzation length [jj, [l6|, ll7|, Q, [l9( have 
yielded reductions of the perturbation result ranging between 7.7 and 8.6 %. In view of 
the simplicity of our analytical treatment leading to a simple transparent picture of the 
band center anomaly in the localization length, we regard the agreement with the numerical 
results of the earlier studies as rather satisfactory. Note also the existence of significant band 
center effects in the coefficients of the quartic terms in (!39l) . as shown by the comparison 
with the corresponding terms in the perturbation result (Eqs. ([8]) and fl22l) ) at a generic 
energy. 

The comparison of the variance ( 138|) with the result ( 1271) for a generic energy reveals the 
existence of two main types of non-perturbative band center anomalies: an enhancement of 
9,1 % of the coefficient of the v'^N term, on the one hand, and the existence of an additional 
weak quadratic term proportional to {Nv'^Y ~ ((~ l^^Af))^) on the other hand. Such a 
quadratic term at 4th order is obtained here for the first time. This term is comparable in 
magnitude with the leading second order term for values of Nv'^ of the order of 10 or larger. 
However, the addition of this term alone does not spoil the single-parameter scaling obtained 
when restricting to the second order in the disorder. Finally it appears that the accurate 
description of the Kappus-Wegner anomaly in the localization length, using the generalized 
improved Born approximation, is crucial for studying the variance since the magnitude of 
the new A^^-term in fl38|) is directly related to this anomaly. 

We close this section with a brief remark about the effect of the non-linear term 

which has been omitted on the r.h.s. of f|T9|) in the expansion of ([3]) for weak disorder. For the 
purpose of a crude estimate we approximate this term by the linearized form a{Rn-i)Rn-i, 
where (Rn-i) is given by ( l30l) . and solve ( fT9i) in the presence of this approximate additional 
term. Thus we replace gn in (|T9l) by 

9'n = 1 - ^I'n^^] , (41) 



1 + ^7. 

at the band center. In this approximation the dominant effect of the non-linear term in the 
moments fllOtilip of — Int^v arises via the denominator 1 — C = 1 — {{g'n)^) in fl35l) which is 
responsible, in particular, for the Kappus-Wegner anomaly in the localization length. From 
the expansion of the parameters in (HTl) for weak disorder it follows that the correction term 
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in this expression leads to a 4th order correction proportional to (z^^)^ in 1 — C, thus leaving 
the dominant second order term of 1 — C unchanged. Our estimatethus shows that the non 
linear term ( l40l) has no effect on the mean and on the variance of — In^Ar to leading order 
in the disorder. 

Finally, we note that a more realistic study of the term fHUl) might be to linearize it in terms 
of the exact solution of the improved Born relation f|T6l) given by 

n . n ^ 

R^ = y e^iHn-m)^^ TT (42) 

m=l ™ p=m+l P 

However, the analytic solution of ( |T9l) in the presence of such a linearized form of ( 140|) is 
clearly very complicated and will not be discussed. 



IV. CONCLUDING REMARKS 



In this paper we have discussed a new approach for studying conductance cumulants for 
weak disorder in the Anderson model in the localized regime. It is based on determining 
complex reflection amplitudes using linear approximations of the exact non linear recursion 
relation between the reflection amplitudes i?„ and Rn-i of disordered chains of lengths na 
and {n — l)a, respectively. These approximate linear relations differ, however, essentially 
from standard weak disorder (perturbation) expansions in that they correctly account for 
the asymptotic {n — > cxd) unitarity property of the averaged reflection coefficient (|i?rip), in 
the localized regime, for any strength of a finite disorder. 

Our analysis relates various anomalous effects of the disorder in the mean and the variance 
of — In in the localized regime to the above asymptotic behavior of the mean reflection co- 
efficient and/or to corresponding behavior of the second moments of the reflection coefficient 
and of the complex reflection amplitude (at the band center), respectively. This includes the 
well-known proportionality of var (— Intjv) to (— Int^v) ~ '^'^N (to leading order) and the 
existence of a leading numerically constant offset in var (— In^Ar) at a generic energy. It also 
includes both the Kappus- Wegner \v\ anomaly and a leading constant offset in (— Int^v) at 
the band center, as well as similar band center anomalies in var (— In^Ar), in particular the 
existence of a new term proportional to A^^ (which we have obtained for the first time). 

An important aspect of our treatment is that it does not rely on the improper use of 
assumptions about phases such as the phase randomization assumption which has frequently 

14 



been invoked in previous work 

The results of Sect. 3 indicate that beyond second order in the disorder and ignoring 
the offset terms, the variance of — Int^r cannot be expressed in terms of the mean alone. 
This rules out single parameter scaling of the distribution of — Intjy even if one assumes 
the higher cumulants to be negligible (in which case the distribution would be lognormal). 
The third, fourth . . . cumulants can, of course, also be studied using our general approach, 
but, as shown in Sect. 1, this requires perturbation expansions to 6th, 8th . . . order in the 
disorder. 

Another application of the analysis of this paper would be the study of conductance 
cumulants in coupled two- and three-chain systems i.e. for few channel quasi-one dimensional 
systems. Localization in such systems has recently been discussed for weak disorder both 
in the case where all the states at the fermi energy belong to conducting bands of the pure 
few-channel system and in the case where, on the contrary, the states for some of the bands 
correspond to imaginary wavenumbers (evanescent states) at the fermi energy 24]. 
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